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Abstract 

(N ' 

I We investigate the thermal quantum and total correlations in the anisotropic XY spin chain in transverse field. While 
.we adopt concurrence and geometric quantum discord to measure quantum correlations, we use measurement-induced 
C I nonlocality and an alternative quantity denned in terms of Wigner-Yanase information to quantify total correlations. 
^> .We show that the ability of these measures to estimate the critical point at finite temperature strongly depend on the 
O anisotropy parameter of the Hamiltonian. We also identify a correlation measure which detects the factorized ground 
1 ^ state in this model. Furthermore, we study the effect of temperature on long-range correlations. 
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r~| '1. Introduction 

The notion of entanglement, on top of being consid- 
£^ ered as the resource of quantum computation, quantum 
cryptography and quantum information processing, has 
^ ,alsn proved to be very useful in analyzing the behavior 

i__iof various condensed matter systems [1, 2]. However, en- 
tanglement is not the only kind of meaningful correlation 

\l present in quantum systems. Quantum discord [3, 4], de- 
fined as the discrepancy between quantum versions of two 
classically equivalent expressions for mutual information, 
has been demonstrated to be a novel resource for quan- 
turn computation [5-7]. Following the path paved by the 

\Q authors of Ref. [3, 4], several new quantifiers of quantum 
correlations and non-locality, that are more general than 
entanglement, have been proposed recently [8-11]. 

Quantum phase transitions (QPTs) are sudden changes 
occurring in the ground states of many-body systems when 

KJ{ one or more of the physical parameters of the system are 
continuously varied at absolute zero temperature [12]. Th- 
C3 ese radical changes, which strongly affect the macroscopic 
properties of the system, are manifestations of quantum 
fluctuations. Despite the fact that reaching absolute zero 
temperature is practically impossible, QPTs might still be 
observed at sufficiently low temperatures, where thermal 
fluctuations are not significant enough to excite the sys- 
tem from its ground state. In recent years, the methods 
of quantum information theory have been widely applied 
to quantum critical systems. In particular, entanglement 
and quantum discord (QD) have been shown to identify 
the critical points (CPs) of QPTs with success in several 
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different critical spin chains, both at zero [13-25] and fi- 
nite temperature [26-28]. It has also been noted that un- 
like pairwise entanglement, which is typically short ranged, 
QD does not vanish even for distant spin pairs [18]. 

An interesting aspect of quantum spin chains in trans- 
verse magnetic field is the occurrence of a non-trivial fac- 
torized ground state [29] . In order to gain a complete un- 
derstanding of these factorized states, the effects of spon- 
taneous symmetry breaking (SSB) should be considered 
[30-32]. In fact, concurrence is known to signal the factor- 
ization point of the anisotropic XY chain corresponding to 
a product ground state [32]. Moreover, it has been demon- 
strated that QD is also able to detect such points, provided 
that either SSB is taken into account or QD is calculated 
for different spin distances [33,34]. In the latter case, the 
factorization point appears via a single intersection of the 
curves of QD. 

In this work, we investigate the pairwise thermal quan- 
tum and total correlations in the one-dimensional anisotro- 
pic spin- 1/2 XY chain in transverse magnetic field. As a 
measure of genuine quantum correlations, we utilize entan- 
glement quantified by concurrence [35,36], and a very re- 
cently proposed observable measure (OMQC) [11], which 
is a simplified version of geometric measure of quantum 
discord (GMQD) [8]. This observable measure has the 
advantage that it does not require a full tomography of 
the system, making it experimentally very accessible. On 
the other hand, in order to quantify non-locality or total 
correlations in a quantum system, we use measurement- 
induced nonlocality (MIN) [9], and an alternative new 
measure defined in terms of Wigner-Yanase skew infor- 
mation (WYSIM) [10]. By comparatively studying the 
thermal quantum and total correlations in the parameter 
space of the Hamiltonian for both first and second near- 
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est neighbor spins, we observe that all of these measures 
are capable of indicating the CP of QPT at absolute zero. 
When the temperature is slightly above the absolute zero, 
in the experimentally accessible region, we analyze the 
ability of correlation measures to correctly estimate the 
CP of the transition. It is also remarkable that, among 
the new measures considered in this work, only WYSIM 
is able to identify the factorization point of the XY spin 
chain even if we disregard the effects of SSB. Lastly, we 
study the long-range correlations of the system and the 
effect of temperature on these correlations. 

This paper is organized as follows. In Section 2, we re- 
view the correlation measures used in our work. In Section 
3, we introduce the one-dimensional XY spin chain, and 
discuss the thermal correlations of the system from differ- 
ent perspectives to estimate the CP of the QPT. Section 
4 includes the summary of our results. 

2. Correlation Measures 

In this section, we review several measures of total and 
quantum correlations to be used in our investigation of 
the one-dimensional XY spin chain. We will commence by 
introducing MIN which encapsulates more general kind of 
correlations than quantum non-locality connected with the 
violation of Bell inequalities [9]. It is defined by (taking 
into account the normalization) 



N{p ab ) 



2 max I 



ab 



U a (p 
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where the maximum is taken over the von Neumann mea- 
surements II a = {11^} that do not change p a locally, mean- 
ing njj/9 a n^ = p a , and ||.|| 2 denotes the square of the 
Hilbcrt-Schmidt norm. MIN aims to capture the non-local 
effect of the measurements on the state p ab by requiring 
that the measurements do not disturb the local state p a . 
It is always possible to represent a general bipartite state 
in Bloch basis as 




Although a closed formula for the most general case of 
bipartite quantum systems is not known, provided that 
we have a two-qubit system (m = n = 2), MIN can be 
analytically evaluated as 

(2{*rt-^#TT<g) ifzVO, 
\2(trTT t -A 3 ) iff=0, 

where TT* is a 3 x 3 dimensional matrix with A3 being 
its minimum eigenvalue, and ||x|| 2 = x 2 with x = 
(x\, X2, xs) f . Due to the symmetries of the considered sys- 
tem in this work, the two-spin reduced density matrix is 
Y-shaped 
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Since the local Bloch vector x is never zero in our investi- 
gation, MIN takes the simple form 



(6) 



Very recently, a new measure of total correlations has 
been proposed in Ref. [10] by making use of the notion of 
Wigner-Yanase skew information 



I{p,X) = -\tv[^p,Xf 



(7) 



which has been first introduced by Wigner and Yanase [37] . 
Here X is an observable (an Hermitian operator) and [., .] 
denotes commutator. For pure states, I(p,X) reduces to 
the variance V(p, X) = trpX 2 — (tipX) 2 . Since the skew 
information I(p,X) depends both on the state p and the 
observable X, Luo introduced an average quantity in order 
to get an intrinsic expression 



Q(p) = J2l(p,X i ), 



(8) 



where {Xi} is a family of observables which constitutes 
an orthonormal basis. Global information content of a 
bipartite quantum system p ab with respect to the local 
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where the matrices {Xi : i = 0, 1, • • • , m 2 — 1} and {Yj : 
j = 0, 1, • • • ,n 2 -l}, satisfying ti(X k Xi) = tr(F fe Y ( ) = S k i, 
define an orthonormal Hermitian operator basis associated 
to the subsystems a and b, respectively. The components 
of the local Bloch vectors x = {x^, y = {yj} and the 
correlation matrix T — tij can be obtained as 



Xi = tip ab {X t ®I b )/^ 

yj =trp ab (I a ®Yj)/^, 
Uj =tvp ab (X i ®Y j ). 



Qa{p ab )=Y, I (P ab ' X *® lb )' 



(9) 



which does not depend on the choice of the orthonormal 
basis {Yi}. Then, the difference between the information 
content of p ab and p a ® p b with respect to the local observ- 
ables of the subsystem a can be adopted as a correlation 
measure for p ab , 



(3) 
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where we add a normalization factor 2/3. Despite the fact 
that the evaluation of most of the measures requires a po- 
tentially complex optimization process, F(p ab ) (referred 
as WYSIM) has the advantage that it can be calculated 
straightforwardly. At this point, we note that quantum 
mutual information (QMI) has been widely used as the 
original measure of total correlations contained in quan- 
tum states. Being based on the von Neumann entropy, 
QMI is a well established measure from the communica- 
tion perspective, while WYSIM is based on the skew infor- 
mation and has a fundamental role in quantum estimation 
theory [10] . Although we do not adopt QMI as a measure 
of total correlations in this work, it has been shown to 
identify CPs in quantum critical systems [14,38-40]. 

GMQD has been introduced to overcome the difficul- 
ties in the evaluation of the original QD [8] . It measures 
the nearest distance between a given state and the set of 
zero-discord states. Mathematically, it is given by 



D G (p ab )=2mm\\p ab 



(11) 



where the minimum is taken over the set of zero-discord 
states. In a recent work, Girolami et al. have obtained an 
interesting analytical formula for the GMQD of an arbi- 
trary two-qubit state [11] 

D G (p ab ) = 2(trS - max^}), (12) 



where S 



;£* + TT* and 
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with {q j = {0,27r,47r} and 9 = arccos{(2trS' 3 -9trS'tr > S' 2 + 
9trS ,3 ) v /2/(3trS' 2 - (trS*) 2 ) 3 }. Furthermore, observing that 
cos ( e+ " i ) reaches its maximum for on = and choosing 



9 to be zero, they have found a very tight lower bound to 
the GMQD, given by 



Q{p ab ) = o ( 2tr5 - V^trS 2 - 2(trS) 2 ) 



(14) 



This quantity (referred as OMQC) can be regarded as a 
meaningful measure of quantum correlations on its own 
and it has the desirable feature that it needs no optimiza- 
tion procedure. Besides being easier to manage than the 
original GMQD, it can be measured by performing seven 
local projections on up to four copies of the state. Thus, 
Q(p) is also very experimentally friendly since one does 
not need to perform a full tomography of the state. 

Lastly, we utilize concurrence to quantify the entan- 
glement content of two-qubit density matrices [35,36]. In 
order to evaluate concurrence, one first needs to calculate 
the time-reversed or spin-flipped density matrix p which is 
given by 

p = (a v ®a y )p*{a y ®a y ). (15) 

Here a v is the Pauli spin operator and p* is obtained from 
p via complex conjugation. Then, concurrence reads 



C(p) = max jo, 



Ai — \ ^2 ~ v A3 — v A 



,}, (16) 



where {A^} are the eigenvalues of the product matrix pp 
in decreasing order. For the simple form of the reduced 
density matrix given in Eq. (5), concurrence reduces to 

C = 2max{0, \pu\ - \p22l IP23I - y/pnpu}- (17) 



3. Correlations in the XY Model 

We start this section with the analysis of the thermal 
quantum and total correlations in the one-dimensional XY 
spin chain in transverse magnetic field. The Hamiltonian 
of the model is given by 



A 



N N 

D(i+«i + (i-7W + j-E«j 

. (18) 

where N is the number of spins, c" (a = x,y,z) is the 
usual Pauli operators for a spin-1/2 at jth site, 7 (0 < 
7 < 1) is the anisotropy parameter and A is the strength 
of the inverse external field. For 7 = the above Hamil- 
tonian corresponds to the XX model. When 7 > it is in 
the Lsing universality class, and reduces to the Ising Hamil- 
tonian in a transverse field for 7 = 1. We are interested 
in the region where the XY model exhibits two phases, a 
ferromagnetic and a paramagnetic phase, which are sepa- 
rated by a second-order QPT at the CP A c = 1. In the 
thermodynamic limit, the XY model can be solved exactly 
via a Jordan- Wigner map followed by a Bogoluibov trans- 
formation. Reduced density matrix of two spins i and j 
depends only on the distance between them, r = \i — j\, 
due to the translational invariance of the system. 




Figure 1: The thermal total correlations quantified by MIN and 
WYSIM as a function of A for 7 = 0.001, 0.5, 1 at kT = (solid line), 
kT = 0.1 (dashed line) and kT = 0.5 (dotted line). The graphs are 
for first nearest neighbors. 



The Hamiltonian is also invariant under parity transfor- 
mation, meaning it exhibits Zi symmetry. Taking these 
properties into account, and neglecting the effects of spon- 
taneous symmetry breaking (which are studied in Ref. [30- 
34]), the two-spin reduced density matrix of the system is 
given by [13] 
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where /o,r is the four-dimensional identity matrix. The 
transverse magnetization is given by [41] 
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where w = v /( 7 Asin0) 2 + (1 + A cos (f>) 2 /2, (3 = l/k b T 
with kb being the Boltzmann constant and T is the ab- 
solute temperature. Two-point correlation functions are 
defined as [42] 
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In Fig. 1, we present our results regarding the thermal 
total correlations quantified by MIN and WYSIM for first 
nearest neighbors as a function of A for kT = 0, 0.1, 0.5 and 
7 = 0.001, 0.5, 1. We note that although MIN and WYSIM 
behave in a similar fashion for 7=1, they show qualita- 
tively different behaviors in the case of 7 = 0.001. Namely, 
WYSIM experiences a more dramatic increase about the 
CP A = 1 than MIN, and reaches to a constant value more 
quickly. Furthermore, it is also important to observe that 
as temperature increases, both of the measures cease to 
exhibit a non-trivial behavior about the CP. 

It has been shown that QPTs can be characterized by 
looking at the two-spin reduced density matrix and its 
derivatives with respect to the tuning parameter driving 
the transition [14,15]. Since correlation measures are di- 
rectly determined from the reduced density matrix, they 
provide information about the CPs and the order of QPTs. 



The CP for a second-order QPT at zero temperature is sig- 
nalled by a divergence or discontinuity in the first deriva- 
tive of the correlation measures. If the first derivative is 
discontinuous, then the divergence of the second derivative 
pinpoints the CP [14-16]. In Fig. 2, we plot the derivatives 
of MIN and WYSIM as a function of A for kT = 0, 0.1, 0.5 
and 7 = 0.001,0.5, 1. We observe that both of the mea- 
sures are capable of spotlighting the CP at kT = for all 
values of 7. It is worth to note that with increasing tem- 
perature, the divergence at CP disappears and the peaks of 
the derivatives start to shift. Therefore, the measures lose 
their significance in determining the CP of the transition. 

We now turn our attention to the analysis of ther- 
mal quantum correlations quantified by OMQC and con- 
currence. In Fig. 3 and Fig. 4, we plot these mea- 
sures and their derivatives with respect to the driving pa- 
rameter A for first nearest neighbors as a function of A 
for kT — 0,0.1,0.5. While concurrence suffers a drastic 
decrease as temperature increases, OMQC still captures 
significant amount of correlation, making it more robust 
against thermal effects. It can also be seen that at kT = 
the CP can be detected by analyzing the non-analyticities 
in the first derivatives of the measures. 

Next, we discuss the question of whether the stud- 
ied correlation measures can signal the emergence of non- 
trivial product ground state in the XY spin chain. Despite 
the fact that the ground state of the model is entangled 
in general, for some special values of 7 and A, the ground 
state becomes completely factorized. In particular, except 
the trivial factorization points A = and A — ¥ 00, there 
also exists a non-trivial factorization line corresponding 




Figure 2: The first derivatives of MIN and WYSIM as a function 
of A for 7 = 0.001,0.5,1 at kT = (solid line), kT = 0.1 (dashed 
line) and kT = 0.5 (dotted line). The graphs are for first nearest 
neighbors. 
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to 7 2 + A~ 2 = 1. Accordingly, as can seen from the be- 
havior of concurrence in Fig. 3 for 7 = 0.5, entanglement 
vanishes at A ~ 1.15, which spotlights the occurrence of a 
product ground state. It is shown in Fig. 2 that, unlike 
OMQC and MIN, WYSIM can signal this factorization 
point through a non-analytical behavior in its derivative. 
For QD to identify this point when the distance between 
the spins is fixed, the effects of SSB must be taken into 
account [33,34,43]. Therefore, it is important to recog- 
nize that the calculation of WYSIM between the spins at 
a fixed distance enables us to detect the product ground 
state even in the absence of SSB. 

Having discussed the behaviors of the thermal total 
and quantum correlations, we now explore the ability of 
these measures to correctly estimate the CP of the QPT at 
finite temperature. Despite the disappearance of the sin- 
gular behavior of MIN, WYSIM, OMQC and concurrence 
with increasing temperature, it might still be possible to 
estimate the CP at finite temperature [28] . For sufficiently 
low temperatures, divergent behaviors of the first deriva- 
tives of correlation measures at T — will be replaced by 
a local maximum or minimum about the CP. Therefore, 
in order to estimate the CP, we search for this extremum 
point. On the other hand, a discontinuous first derivative 
at T = requires us to look for an extremum point in the 
second derivative for T > 0. In Fig. 5, we present the 
results of our analysis regarding the estimation of CP as a 
function of kT for first and second nearest neighbors when 
7 = 0.001,0.5,1. Before starting to compare the ability 
of MIN, WYSIM, OMQC and concurrence to indicate the 
CP in detail, we notice that the success rates of these 
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Figure 3: The thermal quantum correlations quantified by OMQC 
and concurrence as a function of A for 7 = 0.001,0.5,1 at kT = 
(solid line), kT = 0.1 (dashed line) and kT = 0.5 (dotted line). The 
graphs are for first nearest neighbors. 




Figure 4: The first derivatives of OMQC and concurrence as a func- 
tion of A for 7 = 0.001,0.5,1 at kT = (solid line), kT = 0.1 
(dashed line) and kT = 0.5 (dotted line). The graphs are for first 
nearest neighbors. 

measures strongly depend on the anisotropy parameter of 
the Hamiltonian. 

In the case of first nearest neighbors, at 7 = 1, all of 
the correlation measures are able to predict the CP re- 
liably, with concurrence being the most effective among 
them. When 7 = 0.5 MIN turns out to be the worst CP 
estimator. While WYSIM and concurrence points out the 
CP relatively well as compared to MIN, OMQC clearly 
outperforms all others and estimates the CP in a excep- 
tionally accurate way. For 7 = 0.001, MIN and OMQC 
become identical, and they predict the location of the CP 
significantly worse than WYSIM and concurrence. 

For second nearest neighbors, even though we do not 
present the graphs of correlation measures and their deriva- 
tives, the CP has been inspected by performing the same 
analysis as in the first nearest neighbor case. The CPs 
estimated by WYSIM, OMQC and MIN for 7 = 1 deviate 
from the true CP by the same amount but they are still 
acceptable. In the case of 7 = 0.5, both concurrence and 
OMQC estimate the CP very well in contrast to WYSIM 
and MIN. Finally, when 7 = 0.001, while WYSIM and 
concurrence spotlight the CP remarkably well, OMQC and 
MIN perform very poorly. It is also worth to notice that 
concurrence performs even better than the first nearest 
neighbors case for 7 = 0.5 and 7 = 0.001. 

Furthermore, the ability of entanglement of formation 
(EOF) and QD to estimate the CP of the XY spin chain at 
finite temperature has been recently studied by Werlang 
et al. [28]. The performance of the correlations measures 
used in this work as compared to QD and EOF depend on 
the anisotropy parameter of the Hamiltonian and also on 
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Figure 5: The estimated values of the CP as a function of kT for three different values of the anisotropy parameter 7 = 0.001,0.5, 1. The 
CPs in the graphs are estimated by OMQC (denoted by o), WYSIM (denoted by +), MIN (denoted by *) and concurrence (denoted by x). 
Concurrence is not included for 7 = 1 and r = 2, since it vanishes at even very low temperatures. 



the distance between the spin pairs. For instance, in the 
first nearest neighbors case at 7 = 0.5, among the corre- 
lation measures considered here, only OMQC performs as 
well as QD and EOF. On the other hand, for the second 
nearest neighbors at 7 = 0.001, while WYSIM and concur- 
rence turn out to be better CP estimators than QD and 
EOF, MIN and OMQC do not perform as well. We lastly 
note that apart from a limited number of special cases, 
QD still proves to be the most accurate CP estimator for 
the anisotropic XY spin chain. 

Inspired by the methods of Rcf. [21], we now examine 
the long-range behavior of the thermal total and quantum 
correlations for the one-dimensional XY model in trans- 
verse magnetic field. While entanglement vanishes for dis- 
tant spin pairs even in the ordered ferromagnetic phase, 
QD has been shown to remain non-zero [18, 21]. Fig. 
6 demonstrates our results related to the dependence of 
MIN, WYSIM and OMQC on the distance between the 
spin pairs at finite temperature, for A = 0.75, 0.95, 1.05, 1.5 
and 7 = 0.001,1. In case of 7 = 0.001, neither of the 
correlation measures remain significant when the distance 
between the spin pairs is increased. We can also see that 
the decay of the correlations hasten when the tempera- 
ture rises. For the Ising model limit (7 = 1), even though 
MIN, WYSIM and OMQC approach to a finite value in 
the ordered phase for sufficiently low temperatures, ther- 
mal effects wipe out the correlations between distant spin 
pairs after a certain temperature. 

4. Conclusion 

In summary, we have discussed the thermal quantum 
and total correlations in the one-dimensional anisotropic 



XY model in transverse magnetic field from several per- 
spectives. We have quantified the correlations using re- 
cently proposed correlation measures such as WYSIM, MIN 
and OMQC, and a well known entanglement measure con- 
currence. Analyzing these measures in the parameter space 
of the Hamiltonian for first and second nearest neighbors, 
we have found that all of the considered measures are ca- 
pable of indicating the CP of the transition. Although 
the interesting behavior of the measures in the vicinity of 
the CP disappears as the temperature increases, for suffi- 
ciently low temperatures, it is still possible to estimate the 
CP by looking at the derivatives of the correlation mea- 
sures. We have observed that the ability of the measures to 
predict the CP strongly depend on the anisotropy param- 
eter 7. For instance, while OMQC spotlights the CP with 
a remarkably high accuracy at 7 = 0.5 for first nearest 
neighbors, it performs very poorly at 7 = 0.001. On the 
other hand, WYSIM points out the CP reasonably well at 
7 = 0.001 for both first and second neighbors. Moreover, 
we have shown that, among the new measures considered 
in this work, only WYSIM is able to identify the factor- 
ization point of the XY spin chain even if we disregard the 
effects of SSB. These results demonstrate for the first time 
that OMQC and WYSIM are relevant quantities for iden- 
tifying CPs in concrete physical problems. Next, we have 
investigated how WYSIM, MIN and OMQC are affected 
as we increase the distance between spin pairs. We have 
found that the case of 7 = 0.001 is more susceptible to 
both increasing distance of spin pairs and thermal effects. 
It would also be interesting to study these measures for dif- 
ferent quantum critical spin systems such as XXZ chain. 
Finally, we leave the discussion of the effects of sponta- 
neous symmetry breaking on these correlation measures 
as a future work. 
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Figure 6: Long-range behavior of the thermal total and quantum 
correlations for 7 = 0.001 and 7 = 1 at kT = 0.1,0.5. The circles, 
squares, diamonds and triangles correspond to A = 0.75, A = 0.95, 
A = 1.05 and A = 1.5, respectively. 
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